Atom cooling with an atom-optical diode on a ring 
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We propose a method to cool atoms on a ring by combining an atom diode -a laser valve for 
one-way atomic motion which induces robust internal state excitation- and a trap. We demonstrate 
numerically that the atom is efficiently slowed down at each diode crossing, and it is finally trapped 
when its velocity is below the trap threshold. 
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PACS numbers: 03.75.Be, 42.50.Lc 

There is currently much interest in controlling the mo- 
tion of cold atoms for further (deeper) cooling, quantum 
information processing, atom laser generation, metrol- 
ogy, interferometry, and the investigation of fundamen- 
tal physical phenomena. Cold atoms are relatively easy 
to produce and offer with respect to other particles 
many possibilities for coherent external manipulation 
with lasers, magnetic fields, or mechanical interactions. 
They may be trapped in artificial lattices or even in- 
dividually, can be guided in effectively one-dimcnsional 
wires, or adopt interesting collective behavior in Bose- 
Einstein condensates; also, their mutual interactions can 
be changed, or suppressed. All this fiexibility facilitates 
the translation of some of the concepts and applications 
of electronic circuits into the atom-optical realm to im- 
plement atom chips, atom circuits, or quite generally 
"atomtronics" [H. In this context, efficient elementary 
circuit elements playing the role of diodes or transis- 
tors need to be developed. In particular, we have pro- 
posed and studied a laser device acting as a one-way bar- 
rier for atomic motion 0, 0, 0, 0, 0] > ^^'^ similar ideas, 
near experimental verification, have been considered by 
Raizen and coworkers for atom cooling 0, i. (These 
one-way models rely on atom-laser interactions in the in- 
dependent atom regime, but there are also complemen- 
tary proposals making use of interatomic interaction for 
"diodic" one-way transport [9}.) In the following we pro- 
pose to combine, within a ring, the diode and a trap, 
to achieve cooling and trapping with phase-space com- 
pression. Ring-shaped traps for cold atoms have been 
proposed or implemented for matter-wave interferometry 
and highly precise sensors llOll . for studying the stability 
of persistent currents Uj, |l4|, sound waves, solitons and 
vortices in Bose-Einstein condensates , collisions 
for coherent acceleration 15, [l^ . production of highly 
directional output beams 17|, or quantum computation 
[3]. (For further applications see [l^.) The ring traps 
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are implemented by magnetic waveguides [IC 
purely optical dipole forces [2l| , magnetoelectrostatic po- 
tentials |22l|. overlapping of magnetic and optical dipole 
traps [23|, or misalignment of countcrpropagating laser 
beam pairs in a magneto-optical trap [ij]. 
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FIG. 1: (a) Setting of an atom diode in a ring, (b) schematic 
action of the different lasers on the atom levels for the two- 
level atom diode plus quenching, and (c) schematic spatial 
location of the different laser potentials and their effect on 
the moving atom. They are all taken as Gaussian functions: 
flp{x) = tlpH{x,xp,a), Wa{x) — Wa 11(3;, a;^ , f q ) , where 

a = 1, 2, T, Q; ai,2 = cr; and n(a;, xo, a) = exp ^"'g^a ^ ^ • 



In this paper, we assume for simplicity tight lateral 
confinement so that motion along the ring is effectively 
one dimensional or, more exactly, circular with a length 
An atom diode followed by a trap are put on the ring (see 
Fig. [1^). The initial state of ground-state atoms will have 
some velocity and position width, but the anti-clockwise 
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moving atoms will be reflected by the diode -which can 
be crossed only in the direction of the arrow- so all atoms 
will eventually approach the diode clockwise. After each 
crossing of the diode, the atom, which is now excited, 
is forced to emit a spontaneous photon (quenched), and 
ends up at the bottom of the trap; the atom has to loose 
kinetic energy to escape from the trap, so it is slowed 
down at every crossing until being finally trapped when 
its velocity is below the trap threshold. The process is 
reminiscent of Sysiphus cooling (23|, a difference being 
that both the transfer from ground to excited state in 
the diode and the quenched decay are highly controlled, 
robust and efficient processes. 

Before looking at the quantum-mechanical description, 
we will examine a simple classical model to estimate the 
time scales and the cooling efficiency for different recoil 
velocities. In this classical toy model, the diode and the 
trap are reduced to a point at position xd , and the initial 
particle positions and momenta are distributed according 
to Gaussian probability distributions corresponding to 
the initial quantum distributions \'^q{x)\'^ resp. |<i>o(fc)|'^ 
(see below). In each classical trajectory a random recoil 
kick is imparted at the diode clockwise passage, as in the 
quantum jump calculation done below. The trajectory is 
"trapped" (and eliminated from the ensemble) when the 
velocity becomes smaller than the threshold imposed by 
the trap depth; otherwise a fixed amount of kinetic en- 
ergy corresponding to the well depth ^mv^ is subtracted 
and the motion continues. The results for the trapping 
probability are shown in Fig. [2k- Random recoil affects 
the result in two ways: higher recoil velocities accelerate a 
rapid initial trapping, but slightly increase the time nec- 
essary for cooling and trapping the complete ensemble. 
The number of diode crossings before the atom is trapped 
for an initial velocity v > and no recoil is given by the 
smallest integer n fulfilling v — nvx < vt- The time un- 
til this particle has been trapped is given by the time to 
reach the diode the first crossing, io = —{xq —xd)/v plus 
the total time for the n rounds, t„ = / X]j=i ~ j'^T)~^ ■ 
For the parameters of Fig. Oa. and t; = wq we get n = 2 
and +tn ~ 41 ms. 



Now we switch to the quantum mechanical description. 
The scheme of the diode used here can be seen in Fig. 
[TJd and c and it has been explained before [3, 0, 0, lo]. 
In brief, there are three, partially overlapping laser re- 
gions: two of them are state-selective mirrors blocking 
the excited (level 2) and ground (level 1) state on the 
left and right, respectively of a central pumping region 
on resonance with the atomic transition. If the atom 
is traveling from the right and the velocity is not too 
high, it is reflected by the state-selective mirror poten- 
tial Wih/2. If the atom is traveling from the left in the 
ground state then it will be pumped to the second level 
adiabatically (so that the process is robust and indepen- 
dent of velocity in a broad velocity interval) and then 
pass the potential Wifi/2. Note that this setting, see 
Fig. [TJa, can be realized by a detuned STIRAP trans- 
fer [sj] with just two overlapping lasers. To avoid back- 
wards motion after the atom has crossed the diode we 
assume a third level which decays to the ground state 
with a decay rate 73 and we add a quenching laser cou- 
pling levels 2 and 3 with a Rabi frequency fig, see Fig. 
[T|d. a novelty with respect to previous diode models is 
the addition of a ground state well overlapping partially 
with the quenching laser region. The effect of this well is 
twofold: it subtracts kinetic energy from the ground state 
atoms trying to escape from it, and it also traps eventu- 
ally the cooled atoms. The corresponding Hamiltonian 
using |1) = (1, 0, 0)^, |2) = (0, 1, 0)^, and |3) = (0, 0, 1)^, 
where T means "transpose" , may now be written as 



H = 
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where p^; is the momentum operator, and Qp{x) is the 
Rabi frequency for the resonant transition. All potentials 
are chosen as Gaussian functions according to the caption 
of Fig. 1. The "velocity depth" of the trap used is ut := 



Wt 



1.8cm/s, and it corresponds to the trap 



depth used in the classical simulation. 



We examine the time evolution by means of a one-dimensional master equation, 
|p = -J[i/,p]--|{|3)(3|,p}+ + 73 fdu^il + u^) exp(z^z.x) |1) (3|p|3) (1| exp (-z^^x) , (1) 

where Vrec is the recoil velocity and x is the position operator. The initial condition is taken as a pure state 
p{0) = |\E'o)(^'o|, namely a Gaussian wave packet (see the ca ption of Fig. 2) . 



The master equation ([T]) is solved by using the quantum tonian Heff ~ H — i|73|3)(3|. For large 73 (see (26j). 
jump approach [27'|. A basic step is to solve a time- 
dependent Schrodinger equation with an effective Hamil- r> / ^ 

(3|*(0) « -i^^{2\^{t)). (2) 
73 
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FIG. 2: (a) Trapping probability in the classical model; (b) 
trapping probability Pt,x in the quantum model {xmin ~ 
10 fj,m and Xmax ~ 200 /im); and (c) trapping probability Pt,v 
in the quantum model. Thick, green, dotted line: Vrec ~ 
(in the quantum model: averaged over A'^ = 200 trajectories, 
ftp = 4 X 10*/s, H/i = #2 = 4 X W^/s); solid, red line: Vrec = 
3.5 cm/s (in the quantum model: averaged over = 190 tra- 
jectories, f2p = 1 X 10^/s, M = 1^2 = 1 X lO'^/s). Common 
parameters: / = 400 /im (— 200/im < x < 200 /im), m mass 
of Neon, initial wave packet *o(a;) = / dk $o(fc)e*'°'^ with 

= ^1^(1'0'0)^^^P [ - ^OT^ - - ko)(xo - 
^toko) - i^Jo^], ko = fvo, Ak = f A«, xo = -200^111, 
vo = 5 cm/s, Ai; = 4 cm/s, to = 1ms; other parameters in 
the classical model: vt = 1.8cm/s, xd = 80 /im; other pa- 
rameters in the quantum model: Wt = — 10^/s, Wq — 10^/s, 
XW2 — — 90/im, xp = — 40/im, xwi = 10 ^tm, xt — SO fim, 
XQ = 100 /xm, a = 15 /xm, ax = 30 /xm, aq = lO/v'2/xm. 

Therefore, the three-level Schrodinger equation can be 
approximated by a two-level one with the effective Hamil- 
tonian 

n _ Pi,^f Wiix) + Wt{x) Qpix) 
nappro. - 2\ np{x) W2{x) - tWQix) 

where Wq = The second element of the 

approach is the resetting operation at each jump, 
e^pl^im^ux) (3|«'(t)) — > (l|*(t)), where u G [-1,1] 
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FIG. 3: Evolution of the probability densities versus time (a) 
in velocity space, (b) in coordinate space; Vrec = 0, the other 
parameters can be found in the caption of Fig. [21 



is chosen with the probability density |(1 -I- w^), 
all other amplitudes are set to zero, and the wave 
function is normalized. Because of Eq. this 
can also be done in the two-level approximation, 
-^v/WQWexp(^^Mx)(2|^'(^)) ^ then 
the second level is set to zero and the wave function is 
normalized. 

We start by looking at the case with negligible re- 
coil velocity, i.e. with w^ec — 0. We calculate 
the trapping probability in coordinate space, Pt.x = 
J^'""" dx {x\pii\x), and in velocity space, Pt.v = 
'^■'^ (^IPiib)- The results are shown in Fig. [^b/c 
(thick green dotted line) averaging over TV — 200 tra- 
jectories; a numerical error defined by the difference of 
the result between averaging over N and A^/2 trajecto- 
ries is also plotted in Fig. [2]3/c. The parameters used for 
the atom diode with Vrec — result in a range for per- 
fect "diodic" behavior —11 cm/s <v < 11 cm/s (defined 
as in Q). Fig. [3] shows the evolution of the probabil- 
ity density in velocity space, p(w) = (f|(pii + ^22)^), 
and in coordinate space, 'pix) = (a;|(pii -f P22)|2^)- The 
' final probability densities can be seen in more detail in 
Fig. m (thick, dashed, green line) . In this figure we may 
verify the occurrence of cooling and phase space compres- 
sion, namely, a narrower distribution both in coordinate 
and velocity space. The final trapping probabilities are 
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FIG. 4: Initial probability density (dotted blue line); final 
probability density ai t — 400 ms for Vrec = (thick, green, 
dotted line), «,ec ~ 3.5cm/s (solid, red line); other parame- 
ters see Fig. (2] (a) in velocity space; (b) in coordinate space. 

Pj. ^ = 0.9838 ±0.0009 and Pt,„ = 0.9809 ± 0.0040 (with 
the errors calculated as before). 

Let us now examine the case with recoil velocity Uvec — 
3.5cni/s > vt- Because the average value of the recoil 
velocities is still zero, the cooling method will still work 
except for a small fraction of atoms which may acquire 
by successive random recoils velocities above the break- 
up threshold of the diode. (The parameters used for 
the atom diode with Vrec = 3.5cm/s result in a work- 
ing range — 17.5cm/s < v < 22cm/s.) The trapping 
probability versus time shown in Fig. ^jp (solid, red line) 
shows anyway a high final trapping probability. The 
final probability densities are shown in Fig. [4] (solid, 
red line). The main peak is comparable with the main 
peak without recoil, i.e. the velocity width of the main 
peak is smaller than the recoil velocity. We have finally 
Pt^ = 0.9544 ± 0.0003 and Pt,v = 0.9654 ± 0.0017. 

In summary, we have proposed and numerically 
demonstrated a method to cool atoms on a ring, even 
below recoil velocity, after repeated passages across an 
atom diode combined with a ground state trap. 
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